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'^ : Abstract 



The renormalization group has played an important role in the physics of 

>' 
Q> ' the second half of the twentieth century both as a conceptual and a calcula- 



x 



tional tool. In particular it provided the key ideas for the construction of a 



(N 

o 

(^ • qualitative and quantitative theory of the critical point in phase transitions 

o' 

and started a new era in statistical mechanics. Probability theory lies at the 
C ' foundation of this branch of physics and the renormalization group has an 

r^ ■ interesting probabilistic interpretation as it was recognized in the middle sev- 

o' 

O ■ enties. This paper intends to provide a concise introduction to this aspect of 



the theory of phase transitions which clarifies the deep statistical significance 



C^ ■ of critical universality. 
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I. INTRODUCTION 

The renormalization group (RG) is both a way of thinking and a calculational tool 
which acquired its full maturity in connection with the theory of the critical point in phase 
transitions. The basic physical idea of the RG is that when we deal with systems with 
infinitely many degrees of freedom, like thermodynamic systems, there are relatively simple 
relationships between properties at different space scales so that in many cases we are able 
to write down explicit exact or approximate equations which allow us to study asymptotic 
behaviour at very large scales. 



The first systematic use of probabilistic metliods in statistical mechanics was made by 
Khinchin who showed, using the central limit theorem, that the Boltzmann distribution 
of the single- molecule energy in systems of weakly correlated molecules is universal, that 
is independent of the form of the interaction, provided it is of short range. In his well 
known book |]T]] he emphasizes that physicists had not fully appreciated the generality of 
probabilistic methods so that, for example, they provided a new derivation (usually heuristic) 
of the Boltzmann law for every type of interaction. Similar remarks apply to the first 
applications of the RG in statistical mechanics. RG was introduced as a tool to explain 
theoretically the universality phenomena, the scaling laws, discovered experimentally near 
the critical point of a phase transition. In the first period RG calculations used different 
formal devices, mostly borrowed from quantum field theory, which gave good qualitative and 
quantitative results. It was soon realized that a new class of limit theorems in probability 
was involved. This class referred to situations of strongly correlated variables to which the 
central limit theorem does not apply, that is situations opposite to those considered by 
Khinchin. In fact, it was discovered that a critical point can be characterized by deviations 
from the central limit theorem. 

Before providing a description of the content of the present paper we give a short account 
of the development of RG ideas in statistical mechanics. 

It is useful to distinguish two different conceptual approaches. The first use of the RG in 
the study of critical phenomena was based on a Green's function approach to statistical 
mechanics which paralleled quantum field theory. We recall Eq. (1) of 

g{x, {y,}, a) = Z{t, {yi}, a) g{x/t, {y^/t}, aZy^{t, {y,}, a)Z'^{t, {yi}, a)), (1.1) 

where ^ is a dimensionless two-point Green function depending on a momentum variable 
X, a set of physical parameters yi and the intensity of the interaction a (all dimensionless). 
This is an exact generalized scaling relation which in the vicinity of the critical point reduces 
to the phenomenological scaling due to the disappearance of the irrelevant parameters. The 
scaling functions Z and Zy can be expressed in terms of the Green's functions themselves 



via certain normalization conditions. This equation provided a qualitative explanation of 
scaling and, after the introduction of a non integer space dimension d and the use oi e = 4 — d 
as a perturbation parameter 0, became the basis for systematic quantitative calculations 

i, i, i- 

The second approach started with the use on the part of Wilson of a different notion of 

RG 0] that he had already introduced in a different context, the fixed source meson theory 

|§|, [^, with no reference to critical phenomena. This was akin to certain intuitive ideas 

of Kadanoff |T^ about the mechanism of reduction of relevant degrees of freedom near the 



critical point. Kadanoff 's idea was that in the critical regime a thermodynamic system, due 
to the strong correlations among the microscopic variables, behaves as if constituted by rigid 
blocks of arbitrary size. In Wilson's approach in fact the calculation of a statistical sum 
consisted in a progressive elimination of the microscopic degrees of freedom to obtain the 
asymptotic large scale properties of the system. 

Formally the Green's function and the Wilson method were very different and in partic- 
ular the first one implied a true group structure while the second was a semigroup. Both 
gave exactly the same results and the problem arose of clarifying the conceptual structures 
underlying these methods. In fact many people were confused by this situation and some 
thought that the two methods had little connection with each other. Actually the possibility 
of different RG transformations equally effective in the study of critical properties could be 
easily understood using concepts from the theory of dynamical systems. The critical point 
corresponds to a fixed point of these transformations and the quantities of physical interest, 
i.e. the critical indices, are connected with the hyperbolic behaviour in its neighborhood. 



which is preserved if the different transformations are related by a differentiable map [jTT 
Still the multiplicative structure of the Green's function RG and the elimination of degrees 
of freedom typical of Wilson's approach did not appear easy to reconcile. The formal connec- 
tion was clarified in |T^ where it was shown that to any type of RG transformation one can 
associate a multiplicative structure, a cocycle, and the characterizing feature of the Green's 
function RG is that it is defined directly in terms of this structure. In the probabilistic 



setting the multiplicative structure is related to the properties of conditional expectations 
as discussed in [^ and illustrated in the present paper. The relationship between the two 
approaches is an aspect that has not been fully appreciated in the literature and even an 
authoritative recent exposition of the RG history seems to suggest the existence of basic 
conceptual differences [|T3|. In particular Fisher discusses whether equations like ( |1.1D did 



anticipate Wilson and concludes negatively. Different interpretations in the history of sci- 
entific ideas are of course legitimate but after thirty years of applications of RG in critical 
phenomena and other fields this conclusion does not appear justified. In his 1970 paper |^ 
on meson theories Wilson had emphasized analogies of his approach with the Green func- 
tion RG of Gell-Mann and Low even though a detailed comparison was not yet available. 
A balanced presentation of the different RG approaches to critical phenomena can be found 



in the second edition of |jT^ of which, to my knowledge, there is no English translation. 

More difficult was to understand at a deeper level the statistical nature of the critical 
universality. After 25 years I still think that the language of probability provides the clearest 



description of what is involved. The first hint came from a study by Bleher and Sinai |TB 
of Dyson's hierarchical models where they showed that at the critical point the increase 
of fiuctuations required a normalization different from the square root of the number of 
variables in order to obtain a non singular distribution for sums of spin variables (block 
spin). This normalization factor is directly related to the rescaling parameter of the fields 
in the RG. The limit distribution could be either a Gaussian as in the central limit theorem 
(CLT) or a different one which could be calculated approximately. The next step consisted 
in the recognition that new limit theorems for random fields were involved [|l^, jl^, [18 



The random fields appearing in these limit theorems have scaling properties and some 
examples had already appeared in the probabilistic literature. However these examples fl^ 



were not of a kind natural in statistical mechanics. The new challenging problem posed by 
the theory of phase transitions was the case of short range interactions producing at the 
critical point long range correlations whose scaling behaviour cannot be easily guessed from 
the microscopic parameters. A general theory of such limit theorems is still missing and so 



far rigorous progress has been obtained in situations which are not hierarchical but share 
with these the fact that some form of scahng is introduced from the beginning. 

The main part of the present article will review the connections of RG with limit theorems 
as they were understood in the decade 1975-85. The justification for presenting old material 
resides in the fact that these results are scattered in many different publications very often 
with different perspectives. Here an effort is made to present the probabilistic point of view 
in a synthetic and coherent way. Section X will give an idea of more recent work trying to 
extract some general feature from the hard technicalities which characterize it. 

II. A RENORMALIZATION GROUP DERIVATION OF THE CENTRAL LIMIT 

THEOREM 

CLT asserts the following. Let .^i, ^2) • • • ? ^n, • • • be a sequence of independent identically 
distributed (i.i.d.) random variables with finite variance a^ = E(,^j — E(,^j))^, where E means 
expectation with respect to their common distribution. Then 

El (6 - E(6)) n 



an 



1/2 



iV(0,l) (2.1) 



where the convergence is in law and A^(0, 1) is the normal centered distribution of variance 
1. To visualize things consider the random variables ^j as discrete or continuous spins 
associated to the points of a one dimensional lattice Z and introduce the block variables 
C = 2-"/^ Ef i^ and C = 2-"/^ E?:+i i^ Then 

C„+i = -^(C + a. (2.2) 

Therefore we can write the recursive relation for the corresponding distributions 

Pn+i{x) = V2dy Pn{V2x - y)pn{y) = {TZpn){x). (2.3) 



The non linear transformation 7^ is what we call a renormalization transformation. Let us 
find its fixed points, i.e. the solutions of the equation Tip = p. An easy calculation shows 
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that the family of Gaussians 



''°'"<"' = vh''^ <'^^' 



are fixed points. To prove the CLT we have to discuss the conditions under which the 
iteration of TZ converges to a fixed point of variance o"^. The standard analytical way is to 
use the Fourier transform since 7^ is a convolution. In view of the subsequent developments 
here we shall illustrate the mechanism of convergence in the neighborhood of a fixed point 
from the point of view of nonlinear analysis. There are three conservation laws associated 
with TZ: normalization, centering and variance. In formulas 

Pn+iix)dx = / pn{x)dx, (2.5) 



xpn+i{x)dx = / xpn{x)dx, (2.6) 



x'^Pn+i{x)dx = / x^pn{x)dx. (2.7) 

Therefore only distributions with variance o"^ can converge to a Gaussian pG,cr{x)- We fix 
a = 1 and write pc for pg,i- Let us write the initial distribution as a centered deformation 
of the Gaussian with the same variance 

Pr,{x)=PG{x){l + vh{x)) (2.8) 

where 77 is a parameter. The function h{x) must satisfy 

PG{x)h{x)dx = 0, (2.9) 

PG{x)xh{x)dx = 0, (2.10) 

PG{x)x'^h{x)dx = 0. (2.11) 

Suppose now rj small. In linear approximation we have 

{npr,) = Pg{1 + Vi^) + 0{v'), (2.12) 



where L is the hnear operator 

{LK){x) = 27r-i/' f dye-y'h{y + x2-^/^). (2.13) 

The eigenvalues of C are 

Afc = 2i-'=/2 (2.14) 

and the eigenfunctions the Hermite polynomials. The three conditions above on h{x) can 
be read as the vanishing of its projections on the first three Hermite polynomials. 

The mechanism of convergence of the deformed distribution to the normal law is now 
clear in linear approximation: if we develop h in Hermite polynomials only terms with 
k > 2 will appear so that upon iteration of the RG transformation they will contract to zero 
exponentially as the corresponding eigenvalues are < 1. 

To complete the proof one must show that the non linear terms do not alter the conclu- 
sion. This is less elementary and will not be pursued here. 

A terminological remark. The Gaussian is an example of what is called in probability 
theory a stable distribution. These are distributions which are fixed points of convolution 
equations and, with the exception of the Gaussian, have infinite variance. 

III. HIERARCHICAL MODELS 

Suppose now that the ^j are not independent. A case which has played a very important 
role in the development of the RG theory of critical phenomena is that of hierarchical models. 
To keep the notation close to that of the previous section we write the recursion relation 
connecting the distribution at level n to that at level n + 1 

Pn+l{x) = {iZpn){x) = gn{x^){1lPn){x), (3.1) 

where gn{x'^) is a sequence of positive increasing functions and TZ has the same meaning as 
in the previous section. It is clear that such a dependence tends to favor large values of the 
block variable x and therefore values of the .^j's of the same sign. We call this a ferromagnetic 
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dependence. We make the following choice gn{x^) = -L„e^'-'^/^^"^ , where the constant L„ is 
determined by the normalization condition. This type of recursion arises from the following 
Gibbs distribution 

2" 

diJ = Z-'e-^''-^^'--^'"^lldpo{x,), (3.2) 

1 

where Hn has the following hierarchical structure 

/'" xX 

Hn{Xi,... ,X2n) = Hn-l{Xi, . . . .Xan-i) + Hn-l{x2r^-l+l, . . . ,X2n) " c" f ^ ^ 1 , (3.3) 

Hq = and po{x) is a single spin distribution which characterizes the model. The constant c 
satisfies 1 < c < 2. For c < 1 the model is trivial while for c > 2 becomes thermodynamically 
unstable. 

To understand what happens in the case of dependent variables let us consider the 
hierarchical model defined by a Gaussian single spin distribution where the iteration can be 
performed exactly, that is 

Po{x) = {2n)-'/'e-^'/\ (3.4) 

(7^>o)(:r) = {2nalr'/'e-'^'/'-\ (3.5) 

al=(l-2pJ2(^/2A . (3.6) 

We see that the conservation of variance does not hold anymore under the transformation 
TZ and in fact the variance increases at each iteration and is j3 dependent. When n tends to 
infinity the iteration converges to the distribution 

p{x) = (27r)-i/2(l - 2/?c/(2 - c))i/2e-{i-2/3c/(2-c)).V2 (3,7) 

provided P < Per = 1/c— 1/2, that is the CLT holds if the temperature is sufficiently large. 
At /5 = Per the variance of the limit distribution explodes which means that the fluctuations 
increase faster than (9(2"/^). We try a new normalization of the block variable and consider 
^^ ,^j/2"'c^"'/^. The recursion for the distribution of this variable is 

Pn+i{x) = ly" j dypn{2/^'^x - y)pn{y). (3.8) 



We now follow the same pattern as in the previous section: calculate the fixed points and see 
whether they admit a stable manifold or, in probabilistic language, a domain of attraction. 
The fixed points are the solutions of the equation 

p{x) = Le^^' I dyp{2/c^/^x - y)p{y) (3.9) 

with L determined by normalization. A Gaussian solution is easily found 

Pg{x) = {ao/Trf'e-'''^" (3.10) 

with ao = c/3/(2 — c). We shall again discuss stability in linear approximation by considering 
a centered deformation of ( p.lOp Pr,(x) = pg{x){1 + rjh^x)). For small rj 



inp,)=PG(l + vm) + 0(rj'). (3.11) 

where C is the linear operator 

{jCh){x)= I dye-^''^y\h{{x + y)/^/^) + h{{x-y)/c^/^)) (3.12) 



The corresponding eigenvalues for even eigenfunctions are X2k = 2/0^^, with fc = 0, 1, 2, 3, . . . , 
and the eigenfunctions are rescaled Hermite polynomials of even degree. 

Here h{x) must be considered as the effect of many iterations starting from some initial 
distribution which characterizes the model and is therefore dependent on jS. We now see 
that for 2 > c > 2^/^ the eigenvalues Aq and Ai are > 1. The projection of h over the 
constants vanishes due to the normalization condition so that for the iteration to converge 
we have to impose the vanishing of the projection of h on the second Hermite polynomial. 
In view of the previous remark this will select a special value (3cr, the critical temperature 
for the model considered. In conclusion, for 2 > c > 2^/^ and /5 = Per the fixed point ( |3.10| ) 
has a non empty domain of attraction. 

When c < 2^'^ the Gaussian fixed point becomes unstable and we must investigate about 
the existence of other fixed points. Bifurcation theory tells us that most likely there is an 
exchange of stability between two fixed points and we should look for the new one in the 
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direction which has just become unstable. For c < 2^/^ we have A4 > 1 so the instabihty is 
in the direction 1/4, the Hermite polynomial of fourth order. Define e = 2^/^ — c and look, 
for e small, for a solution of (|3.9|) of the form 



PNGix) = J9g(x)(1 - eaH.i^x)) + 0{e') ^ e-''^^)-'l^-'^^)-'l\ (3.13) 

where r* and u* are the fixed point couplings. The analysis of this case is considerably more 
complicated and gives the following results: the linearization of the RG transformation 
around ( |3.13| ) has only one unstable direction so that by requiring the vanishing of an 
appropriate projection along this direction we obtain a non empty domain of attraction for 
some/5cr. 11, [H- 



IV. EIGENVALUES OF THE LINEARIZED RG AND CRITICAL INDICES 

We illustrate the interpretation of the eigenvalues of the linearized RG at a fixed point 
in the context of hierarchical models, which is especially simple. Notations are as in the 
previous section. Consider for definiteness the region 2^/^ < c < 2^/^ so that the Gaussian 
fixed point is unstable, but in its neighborhood there exists a non trivial non Gaussian 



fixed point of the form (|3.13|) with a non empty domain of attraction. Suppose now that 
we start the iteration of the RG from some initial distribution which is close to it but not 
in its domain of attraction. For example we may consider a distribution p(x,/3) of the 



form (|3.13| ) with parameters r, m slightly different from r*,M* which are the values taken 



at the critical temperature /3cr- Application of the RG transformation will eventually drive 
this distribution away from the fixed point due to the presence of an unstable direction. 
However we can "renormalize" our parameters r, u by compensating the instability at each 
iteration and find a sequence of parameters r„, Un such that when n tends to infinity we have 
a sequence of distributions p„(x,/5„) approaching a definite limit. Since we have assumed 
that the parameters r„, «„ are close to the fixed point values, the renormalized parameters 
can be simply expressed in terms of rescalings determined by the eigenvalues of a linear 
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operator analogous to L introduced in the previous section in connection with the Gaussian 
fixed point. This can be seen as follows. Let us write our initial distribution as a deformation 
of (iH) 



p(x,/?) =]9ArG(a;,/3cr)(l + ^/i(a;,/?)). 



(4.1) 



If we develop h in terms of eigenfunctions of the linearized RG at (|3.13|) , the iteration of the 
RG transformation will multiply the projections of h along these eigenfunctions by powers of 
the corresponding eigenvalues. The explosion in the unstable direction can then be controlled 
by rescaling at each step the projection in this direction with a factor proportional to the 
inverse eigenvalue. 

Let us define the susceptibility of a block of 2" spins 

2" 



X„(/?) = ^E 



E^^ 



(4.2) 



This quantity can be easily expressed in terms of the distribution Pn(a;, /3) of the block with 
the critical normalization 2"c~"/^ 

X„(/5) = {2/cr fdxx^pnix,(3). (4.3) 

As n — »• cxD, Xn diverges if /3 = Per- 

To calculate the susceptibility critical index as (3 approaches the critical value we assume 
that the two limits n ^ oo and (3 — >■ (3cr can be interchanged so that they can be calculated 
over subsequences. We want to compute 

logXn(/5n) 



ly^ 



lim 



-n log (c/2) + log / dy y^Pn{y, Pn 



(4.4) 



log I Pn - Per I "^°o log | Pn " Per \ 

If we now choose |/3„ — Pcr\ ~ A"**, where A is the eigenvalue corresponding to the unstable 
direction of the RG linearized at the fixed point, the integral appearing in this formula will 
be almost constant and we obtain 

log(c/2) 



log A 



(4.5) 



12 



Similar calculations can be done for other thermodynamical quantities like the free energy 
or the magnetization. 

We can summarize the situation as follows: given a model defined by an initial distribu- 
tion po, for P < Per we expect the CLT to hold. For /3 = Per by properly normalizing the 
block variables we have new limit theorems where the limit law has a domain of attraction 
which is a non trivial submanifold in the space of probability distributions called the critical 
manifold. If we start from a distribution which is not in the domain of attraction of a given 
fixed point, but not too far from it, it can still be driven to a regular limit by rescaling at 
each step its coefficients in a way dictated by the fixed point. This defines the so called 
scaling limits of the theory associated to a given fixed point. 

For further reading see [^ , p2[ , |2^ . 



V. SELF SIMILAR RANDOM FIELDS 



The notion of self similar random field was introduced informally in |]T6[ and rigorously 



in Wjl and independently in 18 1. It was then developed more systematically in [M and 



25| . The idea was to construct a proper mathematical setting for the notion of RG a la 
Kadanoff- Wilson. This led to a generalization of limit theorems for random fields to the 
situation in which the variables are strongly correlated. 

Let Z'^ be a lattice in (i- dimensional space and j a generic point of Z'^, j = (ji, J2, ■■■,jd) 
with integer coordinates ji. We associate to each site a centered random variable C,j and 
define a new random field 

e; = (7^«,nO. = n-'"/' E ^- (5.1) 



sGVp 



where 

Vp = {s : jkU - n/2 < s^ < Jtu + n/2} (5.2) 

and 1 < q; < 2. The transformation ( |5.1| ) induces a transformation on probability measures 
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according to 

(7^;„/i)(A) = /i'(A) = /i(7^-,U), (5.3) 

where A is a measurable set and 7^* „ has the semigroup property 

7?* 7?* = 7?* fS 4) 

A measure /i will be called self similar if 

K,n/^ = /^ (5.5) 

and the corresponding field will be called a self similar random field. We briefly discuss the 
choice of the parameter a. It is natural to take 1 < a < 2. In fact a = 2 corresponds 
to the law of large numbers so that the block variable ( ^.1|) will tend for large n to zero 
in probability. The case a > 1 means that we are considering random systems which 
fluctuate more than a collection of independent variables and a = 1 corresponds to the 
CLT. Mathematically the lower bound is not natural but it becomes so when we restrict 
ourselves to the consideration of ferromagnetic-like systems. 

A general theory of self similar random fields so far does not exist and presumably is very 
difficult. However Gaussian fields are completely specified by their correlation function and 
self similar Gaussian fields can be constructed explicitly |T^, ||26|. It is easier if we represent 
the correlation function in terms of its Fourier transform 

EteO) = r l[dhp{Xi,...,Xd)e'^^^''^'-'^''. (5.6) 

J —IT 1 

The prescription to construct p in such a way that the corresponding Gaussian field satisfies 
( ^.5|) is as follows. Take a positive homogeneous function /(Ai,...,Arf) with homogeneity 



exponent (i(l + a), that is 

/(cAi,...,cA,) = c<^(i+°)/(Ai,...,Arf) (5.7) 

. Next we construct a periodic function (^(Ai, ..., A^) by taking an average over the lattice 



Z 



d 



^(Ai, ..., A,) = 5^ ^ (5.8) 
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If we take now 

p(Ai, . . . , A,) = n |1 - e^'' l'^(Ai, . . . , A,), (5.9) 



it is not difficult to see that the corresponding Gaussian measure satisfies (|575|). The peri- 
odicity of p insures translational invariance. 

For d = 1 there is only one, apart from a multiplicative constant, homogeneous function 
and one can show that the above construction exhausts all possible Gaussian self similar 
distributions. For c? > 1 it is not known whether a similar conclusion holds. 

From this point on one can follow in the discussion the same pattern as for hierarchical 
models and investigate the stability of the Gaussian fixed points Pq of (|5.5|) . Consider a 
deformation Pg(1 + /i) and the action of T^^n ^^ this distribution. It is easily seen that 

-K^nPch = E(/i|{e;})K,,PG = Hh\{e;})PGm})- (5.10) 

The conditional expectation on the right hand side of ( [5.10D will be called the linearization 
of the RG at the fixed point Pq- To proceed further in the study of the stability we have to 
find its eigenvectors and eigenvalues. These have been calculated by Sinai. The eigenvectors 
are appropriate infinite dimensional generalizations of Hermite polynomials Hk which are 
described in full detail in |2^. They satisfy the eigenvalue equation 



^Hkli^J}) = n['=("/2-i)+il^if,({e;}). (5.11) 

We see immediately that H2 is always unstable while H4 becomes unstable when a crosses 
from below the value 3/2. By introducing the parameter e = a — 3/2, in principle one can 
construct, as in the hierarchical case, a non Gaussian fixed point. The formal construction 
is explained in Sinai's book p^ where one can find also an exhaustive discussion of the 
questions, mostly unsolved, arising in this connection. A different construction of a non 
Gaussian fixed point, for d = 4 has been made recently by Brydges, Dimock and Hurd. This 
will be briefiy discussed in section X. 
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VI. SOME PROPERTIES OF SELF SIMILAR RANDOM FIELDS 

We have already characterized the critical point as a situation of strongly dependent 
random variables, in which the CLT fails. We want to give here a characterization which 
refers to the random field globally. Consider in the product space of the variables ^j the 
cylinder sets, that is the sets of the form 

{eneA,...,e.„GA„}, (6.1) 

with 2i, . . . ,«„ G A, A being an arbitrary finite region in TJ^ and the Ai measurable sets in 
the space of the variables ^j. We denote with Sa the cr-algebra generated by such sets. We 
say that the variables ^j are weakly dependent or that they are a strong mixing random field 
if the following holds. Given two finite regions Ai and A2 separated by a distance 

d(Ai,A2)= min \t - j\, (6.2) 

jeAijGA2 

where |i — j| is for example the Euclidean distance, define 

r(Ai,A2)= sup \^^{AnB)-fI{A)fI{B)\. (6.3) 

AGEA,,BeSA2 

Then r(Ai, A2) -^ when d(Ai, A2) -^ 00. 

Intuitively the strong mixing idea is that one cannot compensate for the weakening of 
the dependence of the variables due to an increase of their space distance, by increasing the 
size of the sets. 

This situation is typical when one has exponential decay of correlations. This has been 
proved for a wide class of random fields including ferromagnetic non critical spin systems 

The situation is entirely different at the critical point where one expects the correlations 
to decay as an inverse power of the distance. In this connection the following result has 



been proved in p8[: a ferromagnetic translational invariant system with pair interactions 



with correlation function 

C{i) = E(eo6) - mo)m^) (6.4) 
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such that 

hm ^^fa^i)<-^<M;^^+i)^(^^ _, (6.5) 

for arbitrary Sk, does not satisfy the strong mixing condition. 

This theorem imphes in particular that a critical 2-dimensional Ising model violates 
strong mixing. Therefore violation of strong mixing seems to provide a reasonable char- 
acterization of the type of strong dependence encountered in critical phenomena. On the 
other end, under very general conditions, if strong mixing holds the one-block distribution 
satisfies the CLT [ ^ . 

An interesting question is whether we can describe the structure of the limit one-block 

distributions that can appear at the critical point beside the Gaussian. It was shown in 

28| , building on previous results by Newman, that for ferromagnetic systems the Fourier 



transform (characteristic function in probabilistic language) of the limit distribution must 
be of the form 

E(e'*«) = e"'*' JJ(1 - tVa|) (6.6) 

j 

with^ . 1/a^ < oo. In the probabilistic literature these distributions are called the D-class 
||36|| . The Gaussian is the only infinitely divisible distribution belonging to this class. 

VII. MULTIPLICATIVE STRUCTURE 

In this section we show that there is a natural multiplicative structure associated with 
transformations on probability distributions like those induced by the RG. This multiplica- 
tive structure is related to the properties of conditional expectations. We use the notations 
of section V. Suppose we wish to evaluate the conditional expectation 

mm}), (7.1) 

where the collection of block variables ^" indexed by j is given a fixed value. Here /i is a 
function of the individual spins C,i- It is an elementary property of conditional expectations 
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that 

E(E(/.|{e7})i{er}) = nh\{^r})- (7-2) 

Let P be the probabihty distribution of the ^j and 7^* ^ the distribution obtained by applying 
the RG transformation, that is the distribution of the block variables ^". By specifying in 
i \l.2\) the distribution with respect to which expectations are taken we can rewrite it as 

E^,„p(Ep(/.|{e;})|{er}) = Ep(/^|{er})- (7-3) 

This is the basic equation of this section and we want to work out its consequences. For this 
purpose we generalize the eigenvalue equation ( ^.11|) to the case in which the probability 



distribution is not a fixed point of the RG. In analogy with the theory of dynamical systems 
we interpret the conditional expectation as a linear transformation from the linear space 
tangent to P to the linear space tangent to 7^^ „-P and we assume that in each of these spaces 

equation pT| 



there is a basis of vectors H[, H, "'" connected by the following generalized eigenvalue 



^piH^m) = A.(n,P)<-"'^({e;}). (7.4) 

Equation (|7.3|) implies that the A's must satisfy the relationship 

A,(m, K,n^)Afc(ri, P) = Xkimn, P). (7.5) 

From ( |7.4| ) and (|7.5|) we find that the A^ are given by the following expectations 

A,(n, p) = mT'-^'^mm^iiQ)), (7.6) 

where H^ are dual to H^ according to the orthogonality relation J H[H^dP = 5k j- The 
Afc are therefore special correlation functions. The similarity between equation ( |7.5| ) and 
( p..l| ) is then obvious. The Green's function RG corresponds to a very simple transformation 
on the probability distribution such that its form is unchanged and only the values of its 
parameters are modified. 



VIII. RG AND EFFECTIVE POTENTIALS 

In this section we want to illustrate a connection between RG and the theory of large 
deviations |3^. By large deviations we mean fluctuations with respect to the law of large 
numbers, e.g., fluctuations of the magnetization in a large but finite volume. In view of the 
connection of RG with limit theorems our discussion will parallel, actually generalize, some 
well known facts in the theory of sums of independent random variables. This will lead to a 
probabilistic interpretation of a widely used concept in physics, the effective potential, and 
will clarify its relationship with the effective Hamiltonian in RG theory. We continue with 
our model system of continuous spins C,i indexed by the sites of a lattice in d dimensions and 
try to estimate the probability that the magnetization in a volume A be larger than zero at 
some temperature above criticality. From the exponential Chebysheff inequality we have for 
^> 0,x> 

P ( J]6/|A| >x]< e-I^I^^E(e^S»^A«0 < e-'^'^d^'-^), (8.1) 

where 

r(|A|,x) =sup(fe-^logE(e^^'eA€»)) (8.2) 

6»>0 |A| 

is the Legendre transform of 4t logE(e^^»6A'^i). With some more work one can estabhsh 

also a lower bound 

P (j2^^/\M >^)> e-l^l(r(l^l'-)+"(l^l)+^), (8.3) 

VieA / 

with a — > for A ^ oo, and 6 > arbitrarily small. We then conclude 

-\im^^\ogp(J2^^/\M>A =2imr(|A|,x) = K//(x), (8.4) 

where Veff{x) is known in the physical literature as the effective potential. An important 
remark. While r(|A|, x), being the Legendre transform of a convex function, is always convex 
for any A, this is not the case with — -4^ logP(^^g^^j/|A| > x) = V"(|A|,x) for finite A and 
one has to be careful in interpreting for example results from numerical simulations. 
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To understand the connection with the RG it is convenient to consider first the case 
of independent random variables, that is the situation considered in section II. A classical 
problem in limit theorems for independent variables is the estimate of the corrections to the 
CLT when the argument of the limit distribution increases with n. A well known result in 
this domain is the following [^: suppose we want to estimate Pi^iii/n^^"^ > x) when 



X = 0(71^^"^). Then 

P (j2^^/^'^' > ^V e-"S^r^(^"~''')'= (8.5) 



1 



for n ^ 00 and liiOn^co x^'^^n^^'^^^^^'^ = 0. The function T{z) = Y^'^T^z'^ is the Legendre 
transform of logE(e^^*). The sign ^ has to be understood as logarithmic dominance. If 
X = 0{n^^^), the whole function T contributes and we are back to the large deviation 
estimate at the beginning of this section. 

We expect a result like (|875|) to hold for the one-block distribution in the case of dependent 
variables as in statistical mechanics away from the critical point. We then see that the 
coefficients of an expansion of r(|A|, x) in powers of x determine the corrections to the CLT 
for the one-block distribution. More interesting is the situation at the critical point. Suppose 
first that the one block limit distribution is Gaussian but the normalization is anomalous 
as it is the case in hierarchical models for a range of values of the parameter c. Instead of 
(fj.5D we expect an estimate of the form 



P (y,^^/\M' >^]^ e-I^IS^^^d^D^^'/l^l"")', (8.6) 

with p > 1/2 and lim|A|->ooa^*'*'^/|A| ^ '^ = 0. We see that for the quadratic term to 
survive the coefficient r2 must vanish when |A| ^ 00 as |A| ^ ^. If the one-block limit 
distribution is not Gaussian we can establish a general relationship between its logarithm 
and the effective potential. Let us write — logP = Vrg where P is the limit distribution. 
We now rewrite the large deviation estimate in the following way 

P ( J^e^/IAr > x\A\^-A ^ e-l^l^d^l'"). (8.7) 
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Scale X ^> xl\h\ ^ . Then we obtai 



am 



VnG{x)= \m^ |A|r(|A|,x/|A|^-0. (8.8) 

Therefore r(|A|,a;) determines in different hmits either V"e// or Vrq- The discussion in the 
present section can be made rigorous in the case of hierarchical models. 

IX. COEXISTENCE OF PHASES IN HIERARCHICAL MODELS 

In the case of hierarchical models the RG recursion relation for the one-block probability 
distribution can be easily rewritten as a recursion for the quantity V^(|A|,x) introduced in 
the previous section which coincides with the effective potential in the limit |A| -^ oo. In fact 
if we normalize the block-spin with its volume, that is consider the mean magnetization, a 
simple calculation gives the following iteration for the corresponding probability distribution 
7r„(a;): 

7r„(x) = L^e^''"'' / rfxV„_i(2a; — x')7r„_i(x'). (9.1) 

Taking the logarithm and dividing by the number of spins 2", we obtain 

Vn{x) = -1/2'' login - /3(c/2)V - 1/2" log /" c^^'e-2""'(^-i (2— ')+v;.-i (-')). (9.2) 

To illustrate the difference between ( p.2| ) and ( |3.8| ) let us consider again the simple case in 
which the model is defined by a Gaussian single spin distribution. The iteration of ( p.2| ) 
gives for small (3 



Vn{x) = 1/2 



l-/3^(c/2)^' 



x' + Un (9.3) 



where Ur, tends to zero when n -^ 00. In this limit then 



■'n 



Veff = l/2[l-2cP/{2-c)y. (9.4) 

The critical temperature is defined as the value (3cr for which the coefficient of x^ vanishes 
and coincides with that found in section III. On the other hand in that section it was the 
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only temperature for which the recursion (|3.8|) converges to the Gaussian fixed point, i.e. 
the only temperature for which the following difference between two diverging expressions 
converges 

n 

2(3j2iVcf-i2/cr^\ (9.5) 



We want to apply now (|9.2|) to the study of the magnetization in the phase coexistence 
region for a general hierarchical model [^, ^^. The problem we want to discuss is the 



following. In the hierarchical case at level n we have blocks, containing each 2" ^ spins, 
interacting in pairs through the Hamiltonian 



c"(C-i + C-i)74 = c-iCnY (9.6) 

where (n-iXn-iXn are mean magnetizations. Suppose now that Cn is assigned the value 
aM, < a < 1, M being the spontaneous magnetization corresponding to the temperature 
/3. We want to calculate the conditional distribution of (n-i given (n, for large n. The 
remarkable result is that one of the quantities Cn~i o^ Cn-i with probability close to 1 is 
equal to the full magnetization M. 

To compute the desired distribution we have to estimate 7r„(x) or, what is the same. 



Vn{x) for large n. From ( |9.2| ) we expect asymptotically 

Vn{x) = Veffix) + {c/2fY{x) + ... (9.7) 

Since in the phase coexistence region Veff{x) is fiat, i.e. constant, the whole x dependence 
is given by Y{x). In order to compute this function we perform a subtraction and consider 
Vn{x) — Vn{xo) choosing xq in the fiat region of Veff{x). From ( p^) it is easily seen that the 
quantity 

A„(x) = (2/c)"(K(a:) - Vn{xo)) (9.8) 

satisfies a recursion of the form 

A„(.) ^ -C-" log^„ - p.' - C-" log / ...e— (--.<-'.«-.^-'.). (9.9) 
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where An is determined by the condition A„(xo) = 0. Let us choose xq = M. By symmetry 
A„(±M) = 0. For < X < M and large n the main contribution to the integral on the 



right hand side of ( |9.9| ) comes from the region x ± x' ~ M, while for —M < x < from the 



region x ± x' ~ —M. We can write therefore the approximate recursion equations 

A„(x) = p{M^ - x^) + c-^A„_i(2x T M), (9.10) 

where the =F in the second term on the right corresponds to < x < M or —M < x < 0. This 



type of equations has been rigorously studied by Bleher |^ and the asymptotic solutions 
show a complicated fractal structure. 

The conditional probability of interest to us is 

f ^3,/g-c"-i(A„_i(x')+A„_i(2Q:M-x')) 

P(|C„-, - Af I < eM|C„ = aM) = "'~p^,^_,.-,,,„_,^,,,„_,^„_,,, (9.11) 

Joo 

Since the main contribution to the integral in the denominator comes from the same region 
appearing in the numerator, our conditional probability is for sufficiently large n as close as 
we want to 1. 

X. WEAK PERTURBATIONS OF GAUSSIAN MEASURES: 
A NON GAUSSIAN FIXED POINT 

Starting at the end of the seventies the RG has become a very important and effective 
tool for proving rigorous results in statistical mechanics and Euclidean quantum field theory. 
An impressive amount of work has been done and it is not possible to give even a schematic 



account of it |^ . Many different versions of the RG idea have been used, each author or 
group of authors following his own linguistic propensities. Probability theory is always in 
the background and we want to try to recover some conceptual feature common to all of 
them. As in the previous part of this review limit theorems will be a relevant reference. 
However the limit theorems to be considered are of a different kind, they are those which in 
probability are called non classical and are related to the following problem. 
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Given a probability distribution P and an integer n, can one consider it as resulting 
from the composition (convolution) of n distributions Pk, k = 1,2, . . . ,n7 In other words 
can one consider the random variable described by P as the sum of n independent random 
variables? In formulas 

P = Pl^P2*...^Pn, (10.1) 

where -k means convolution. It is clear, for example, that a Gaussian distribution of variance 
(T^ can be thought as the composition of any number n of Gaussians with variances af 
provided Yli (^f = o^^- The problem arises naturally of investigating under what conditions 
convolutions like the right hand side of (|10.1|) converge to a regular distribution as n — i> cxd. 



The right hand side of (|10.1|) can be considered as a recurrence relation 

Pn+l=Pn*Pn, (10.2) 



where Pn = Pi-k P2'k . . .*P„_i. Comparing ( |2.3| ) with ( |10.1| ) we see that while in the case of 
limit theorems for independent identically distributed random variables we have a natural 
fixed point problem, this is not in general the case for non identically distributed variables. 
As we shall see below, the RG approach to Euclidean field theory and the statistical mechan- 
ics of the critical point has led to formulations which have analogies with these problems. 
In fact, infinite dimensional equations structurally similar to (|10.2|) are constructed which 
can be transformed into equations admitting fixed points after a rescaling. 

In the following exposition we shall follow the recent article by Brydges, Dimock and 
Hurd [^]. The goal of these authors is the construction of a quantum field theory in R'' with 
non trivial scaling behaviour at long distances, that is in the infrared region, determined by 
a non Gaussian fixed point of an appropriate RG transformation. The starting point is a 0^ 
theory in finite volume regularized at small distances to eliminate ultraviolet singularities. 
This model is believed to have a non Gaussian fixed point in 4 — e dimensions and to simulate 
such a situation in 4 dimensions the authors introduce a special covariance for the Gaussian 
part of the measure. The first step consists in the construction of a covariance v{x — y) 

24 



which behaves at large distances hke (—A) ^ "^Z^. This means that it scales like |x| ^^"^ for 
large |a;|. Their choice is 

/oo 
c/aa^/2-2g-|.-y|V4a_ (10.3) 

Such a covariance can be decomposed in the following way 

v{x -y) = Yl L-^'~'^^C{L-^{x - y)), (10.4) 

j=0 

where L > 1 is a scaling factor and 

C{x) = / rfaa^/2-2g-|.|V4a_ (10.5) 



Each term in the expansion can be interpreted as the covariance of a rescaled field 

(f)L-,{x) = L-(2-^)i/20(^-Ja;) (10.6) 

which has reduced fluctuations and varies over larger distances. The aim is to study the 



measure 

dfiA = Z-^e-^^^'f'Ufi^ (10.7) 

where 

Va(0) = A /:0^:, + C [ :{d<Pf:^ + fx / : 0^ :, (10.8) 

when A tends to R^. The double dots indicate the Wick polynomials with respect to the 
covariance v. 

Take for A a large cube of side L^ so that the measure is well defined. We want to 
calculate 

(/i^^e-^)(0) = (/i^^^/i^^_^^....^/i(5^vce^^)(0) (10.9) 

having used the above decomposition of the covariance with 

Cj{x) = L-^^-'^'C{L-^{x - y)) (10.10) 
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Actually in finite volume we should specify some boundary conditions but we shall ignore 
this aspect. Next by defining 

we find the recursion relation 

4+i(0) = (/x^^^Z,)(0). (10.12) 

In this way the calculation is performed by successive integrations over variables which 
exhibit decreasing fluctuations. This is not yet our RG equation because as j -^ oo, ^(j. 
becomes a singular distribution and we do not obtain a fixed point equation. However by 
introducing the rescaled fields 0L-j(a;) = L^'^^~^)''/^0(L~-'x) and the rescaled Z's 

Z,(0) = Z,(0i-.) (10.13) 

the recursion becomes 

Z,+i(0) = (/i^^Z,)(0^-i) (10.14) 

with initial condition Zq = e~^ . We emphasize that the last step is possible due to the 
special structure of the measures fi^,. It is now meaningful to look for the fixed points of 
( |10.14| ). Brydges, Dimock and Hurd have proved that in c? = 4 for e small there exists a 
non Gaussian fixed point of ( |10.14D characterized by a value A(e, L) of the coupling A and 



that for certain values /i(A), ({X) the iteration of ( |10.14D with initial condition e~^ converges 
to this fixed point. Technically the proof is very complicated and its description is beyond 
the aims of this review. A very good exposition with some simplifications of the techniques 
employed can be found also in 



XI. CONCLUDING REMARKS 

The question we want to consider is the following: which are the benefits for our un- 
derstanding of critical phenomena and more generally of statistical physics deriving from 
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the use of probabilistic language? Feyninan thought that it is worth to spend one's time 
formulating a theory in every physical and mathematical way possible. In our case there is 
an intuition associated with probabilistic reasoning that is foreign to the usual formalisms 
of statistical mechanics based on correlation functions and equations connecting them. 

Apart from this general remark we must consider that the rigorous results obtained so far 
in RG theory have been strongly influenced by the probabilistic language as this appears the 
most natural for the mathematical study of statistical mechanics and Euclidean field theory 
when a functional integral approach is used. New technical ideas however are needed to deal 
with concrete problems like calculating the critical indices of the 3-dimensional Ising model 
or establishing in a conclusive way whether the field theory 0^ is ultraviolet non trivial. 

The formal apparatus of RG has been easily extended to the analysis of fermionic sys- 
tems when these are described by a Grassmaniann functional integral ||39|, that is by the 



analog of a Gibbs distribution over anticommuting variables. In this case the convergence of 
perturbation theory plays a major role on the way to rigorous results. Recently, it has been 
possible to give a true probabilistic expression to general Grassmaniann integrals in terms 
of discrete jump processes (Poisson processes) |^, [^ so that classical probability may 



become a main tool also in the study of fermionic systems especially in view of developing 
non perturbative methods. For an early example of connection between anticommutative 



calculus and probability see [42 



In a wider perspective one may remark that the theory of Gibbs distributions is becom- 
ing instrumental also in various sectors of mathematical statistics, for example in image 
reconstruction, and critical situations appear also in this domain. Transfer of ideas from 
statistical mechanics to stochastic analysis is currently an ongoing process which shows the 
relevance of a language capable of unifying different areas of research. Probability theory 
for a long time has not been included among the basic mathematical tools of a physics cur- 
riculum but the situation is slowly changing and hopefully this will help cross fertilization 
among different disciplines. 
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